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In this talk I will present some of the main difficulties we encounter in studying the large scale 
behavior of disordered systems. This presentation will be done using a field theory language. 
The difficulties in applying the standard renormalization group approach are due to the presence 
of strong non-perturbative effects that we do not know how to master. These difficulties are 
particular acute in the case of ferromagnets in random field, localized electrons, growth models 
and spin glasses. I will review the situation presenting at the end a promising approach based on 
fat diagrams. 
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1. Introduction 

Many disordered systems are characterized by the presence of quenched disorder in the equa- 
tion of motion (e.g. in the Hamiltonian). In other systems, like structural glasses, the Hamiltonian 
is translational invariant and the disorder is only present in the equilibrium configurations [jl|]. Al- 
though the two classes of systems have many points in common, I will consider only systems 
belonging to the first class. In disordered systems there are new phenomena that are not present 
for ordered systems, i.e. disordered systems display a new collective behavior. Sometimes a very 
strong slowing down of the dynamics is present (more than fifteen orders of magnitude): in this 
case the system forms a glass. In these thirty years there have been many progresses and break- 
throughs, however in some cases there are fundamental questions that are still not answered, in 
spite of the fact that some of the questions have formulated twenty years ago (and may be more). 

There are many random systems whose critical properties are not understood. Technically our 
goal is to study the large-scale behavior of random systems using the renormalization group. This 
is a very interesting problem, since random systems have peculiar features (e.g. non-perturbative 
phenomena like Griffiths singularities) and often conventional approaches do not work. There are 
many random systems where the renormalization group approach is not developed enough. 

Disordered systems have also many applications beyond physics. In many case we are inter- 
ested to study the collective behavior of a large heterogeneous system of interacting agents. There 
are many examples in the literature, e.g. morphology of river basins, optimization algorithms, neu- 
ral networks, structure of networks in quite different contests (internet, food, finance, proteins), 
species dynamics, evolution. . . 

In this talk I will concentrate my attention on the unsolved problems. I have tried to construct 
this presentation in a way that Hector would have appreciated; I hope to have reached my aim 
by using as far as possible the language of field theory, that is very compact and Hector mastered 
in a magnificent way. I will be very fast on the motivations for studying the models and on the 
phenomenological consequences of the theories because I want to stress the root of the difficulties. 

This paper is organized as follows. After the introduction in section || I will present some 
general considerations on the renormalization group; in the next section I will show which are the 
main difficulties in applying the renormalization group to some disordered systems. In section 
4 I will discuss some open questions for ferromagnets in a random field; in the next section I 
will discuss the properties of the Schrodinger equation in a random potential. Section 6 will be 
dedicated to spin glasses, in the next section I will present a very interesting new approach, i.e. the 
topological expansion for lattice field theories, where one expands the correlations function using 
fat diagrams. Finally in the last section I will present some brief conclusions. 

2. The renormalization group 

Our aim is to understand the collective behavior of disordered systems, when long-range cor- 
relations are present. This may happens or at the critical point, or in a phase where a continuous 
symmetry is broken and therefore Goldstone Boson type excitations are present. Sometimes we 
are in the strange situation where the symmetry that is spontaneously broken is induced by the 
formalism we use (e.g. replicas). In all these systems we want to get insight, i.e. a quantitative and 
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qualitative understanding. In this respect the study of universal quantities (i.e. those quantities that 
do not depend on the detailed form of the Hamiltonian) play an important role. 

As an example let us consider an homogenous Heisenberg ferromagnet (neglecting long-range 
dipole interactions). For a given material, near the Curie point that separates the high temperature 
phase from the low temperature phase (where spontaneous magnetization is present), the magnetic 
susceptibility x(T) behaves as 

X(T)n —^—- . (2.1) 

The exponent /does not depends on the material, but only the dimensions of the space (e.g. D = 3) 
and it is an universal quantity. On the contrary A and T c do depend on the material. 

The renormalization group is one of the most powerful tools to study the structure of the phases 
of a physical system. Using the renormalization group approach one can obtain both qualitative 
and quantitative results In this way we can construct a general theory of second order phase 
transitions and we can classify the transitions into universality classes. The values of the critical 
exponents can be of two different types: a) they are equal to those of mean field theory as it 
happens above the upper critical dimension (Du); b) they differ from those of mean field theory 
(this happens below the upper critical dimension Dy). It is possible to find out when the phase 
transition disappears together with the low temperature phases. The low temperature phase is 
absent in dimensions below the lower critical dimension (Dl). The exponents are usually a non- 
trivial function of the dimensions between Dl and Dy. One of the great lessons of the seventies 
was that the study of the phase diagram of a system in the temperature-dimension plane is very 
useful and it is a source of inspiration for understanding how a system behave in our three (space) 
dimensional world. The study of the behavior of models in dimensions different from 3 has been 
the starting point of crucial computations in three dimensions. 

The successful quantitative use of the renormalization group for ferromagnetic phase transi- 
tion in 3 dimensions, started by the recognition the 4 is the upper critical dimension (i.e. the <p 4 
theory has a trivial infrared stable fixed point for positive coupling constant in 4 dimension). This 
observation allowed a quantitative control of the infrared stable fixed point (an consequently of the 
critical exponents, i.e. the anomalous dimensions of the operators), leading to the celebrated 4 — e 
expansion for the critical exponents. Later on a similar e expansion was constructed near the lower 
critical dimension: for example in the case of a continuous symmetry breaking of a conventional 
symmetry (e.g. 0(N)) one can derive a 2 + e expansion. 

Let us consider the case of a ferromagnetic Heisenberg model, that may be represented by a 
conventional 4 theory invariant under the internal 0(3) group (this is essentially the Goldstone 
model). We have already seen that the upper critical dimension is 4 and the lower critical dimension 
is 2. In dimensions greater than 4 the theory exists in presence of a cutoff. In this situation infrared 
divergences are absent and one finds the there are no anomaly dimensions. The critical exponent y 
takes its naive mean field value, i.e. 1. In dimensions less or equal to 2 there are no transitions. 

How to compute the value of y at dimensions 3? We can use or a 4 — e expansion or a 2 + e 
expansion (or also an 1 /N expansion at D = 3 for a theory invariant under the 0(N) group). The 
best method is to use a fixed dimensional approach: in this case seven non-trivial terms are known 
(seven loops). In this way one obtains results with and errors of a few parts over 10 . 
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In many systems using the renormalization group one can obtain quantitative prediction on the 
exponents, on the critical equation of state and on other dimensionless quantities. One can reach 
these goals by using an e expansion at dimensions not too far away from the upper or the lower 
critical dimension. If for a given system we are not yet able to compute the values of the upper and 
of the lower critical dimension we miss something important in our understanding; unfortunately, 
as we shall see later, in many random systems we are exactly in this sad situation. 



3. Difficulties with random systems 

The behavior of strongly disordered systems is frequently very different from the behavior 
of pure, homogenous systems: electron localization is a very good example of this fact. Some- 
times the most dramatic effects show up experimentally in a time domain where many different 
phenomena are present: very slow relaxation and aging [Q], memory and oblivion [Q] and gener- 
alizations of the usual fluctuation dissipation relations [|5j]. These dynamical effects have often a 
static counterpart. For example in spin glasses ^ the onset of very slow relaxation and the conse- 
quent divergence of the correlation time is associated to the divergence of a static quantity, i.e. the 
non-linear susceptibility |7|] ; in spin glasses the form of generalized fluctuation dissipation relations 
in the low temperature phase is related to an other static quantity, i.e. the probability distribution 
of the overlap In many cases the starting point of our theoretical approach is the homogenous 
situation where disorder is absent and a perturbation expansion is done in powers of the strength 
of the disorder. It should be not surprising that we have serious difficulties in understanding the 
non-perturbative effects produced by the disorder. 

Sometimes the replica theory is used to average over the disorder. Replica theory [||] is a 
wonderful tool and it allows a very compact reformulation of the problems: often it works in 
a magnificent way and also it suggests new quantities to measure. However sometimes replica 
theory puts the dirty under the carpet and cleaning under the carpet is not an easy job; a strong 
effort must be done to decode replica theory and to give a physical meaning to the results of the 
replica approach. The study of non-perturbative effects within the replica method has a lot of 
subtleties that are only partially understood. 

Due to these difficulties there are many physical and experimental questions on the static 
behavior of random systems that do not have a well established theoretical answer. These questions 
are waiting for an answer from at least twenty years, without getting an unambiguous response. 
Why are we stuck in spite of all progresses that we have done in these recent years? The problems 
arise from our inability to use the renormalization group for random systems in presence of non- 
perturbative phenomena. In many cases the ignorance of the upper (or lower) critical dimension 
forbids us to follow the standard path; this is a clear sign of lack of understanding: indeed in most 
of the cases this is due to the absence of a manageable field theoretical representation that can be 
used as starting point. This is what happens for many disordered systems. I will mention some of 
the most important cases. 



FeiTomagnetic Ising systems in a random field. Here we know both the upper critical dimen- 
sion (Du = 6) and the lower critical dimension (Dl = 2). However dimensional reduction 
[10] (based on supersymmetry), that connects the critical exponents for such a system in 
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dimension D to those of a standard ferromagnet in dimension D — 2, is not valid at the non- 
perturbative level for some known reasons JIo|]. Unfortunately the breaking of dimensional 
reduction is a non-perturbative effect and there is no consensus on how to modify the 6 — £ 
expansion in order to include this effect. Moreover there are some qualitative features of the 
behavior of the model in three dimensions that are not understood (e.g. the existence of a 
field with dimensions very near to zero). 

The localization transition. Here the lower critical dimension is well known {Di = 2), but 
the upper critical dimension is not known. In a field theory approach (where one starts 
from a theory with negative coupling constant) localization is due to non-perturbative 
phenomena and the localization transition is related to the formation of a zero mass bound 



state [11, 12, 13]. Different values of the upper critical dimension have been proposed, but 
no consensus is reached: there are very reasonable indications that the critical exponents are 
non trivial for any dimensions, but the problem should be further analyzed. When electrons 
interact strongly we have Mott localization and we face the same problems. New phenomena 
may appear as an effect of the Mott localization. There is a good experimental evidence that 
in many materials, mainly thin films, slow phenomena set in. May be the localized electrons 



form a kind of glass, the so-called electron glass [14, 15]. Our understanding is limited. 



Spin glasses at zero magnetic field. Here we know the upper critical dimension {Dy = 6) and 
we control the critical exponents in the 6 — £ expansion, but we do not have a good control 
of the infrared behavior in the low temperature phase. There are long-range correlations (i.e. 
massless particles) that are due to the existence of Goldstone modes of a symmetry that is 
spontaneously broken in the low temperature phase. Unfortunately the symmetry group is 
in some sense an infinite dimensional group and already at zero-loops the situation is quite 
complex: the propagators have both poles at k 2 = and at k 2 = m 2 , where the index i goes 



up to infinity and the m ( 's accumulate toward zero [jl6p. We are not able do construct the 
appropriate non-linear sigma-model that should be used near the lower critical dimensions. 
The Ward identities of the symmetry have been identified and they have been shown to be 
the responsible of many important cancellations, but we are still far from writing the final 



effective theory [16]. Usually the non-linear sigma model is obtained from the linear model 
by sending to infinity the mass of the transverse modes in such a way that only the massless 
one (the Goldstone modes) survive. This requires a clear cut-separation among massive and 
massless modes and this is a nightmare when there is an infinite number of particles whose 
mass accumulates to zero. Ironically there is an analytic computation [|l8|] of the lower 
critical dimension that gives the extravagant value Di = 2.5. This value is confirmed by 



numerical simulations Q19| ] and it is based on a non-perturbative computation of the domain 
wall energy: unfortunately we are unable to transport the basic ideas into the perturbative 
field theory formalism based on propagators. 

Spin glasses at non-zero magnetic field. Here the situation is worse: we lack both the lower 
and the upper critical dimensions. At the naive upper critical dimensions (i.e. Dy = 6) 
the trivial fixed is not infrared stable and no infrared stable fixed point is known [^]. The 
infrared stable fixed point could be a strong coupling one, like in four-dimensional QCD, or 
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a zero temperature fixed point; the situation calls for further investigation: many different 
alternative scenarios are possible. 



Surface roughening in the framework or the KPZ equation [ ]21[ ] for surface growth in D + 1 
dimension. (The model can also be mapped on the statistical properties of an equilibrium 
directed polymer in a random medium.) Here the lower critical dimension is known (i.e. 
Dl = 0). The problems come with the upper critical dimension. From time to time arguments 



(albeit not completely convincing) are presented [ ]22| ] that hint at D\j = 4, but this result seems 
to be at variance with numerical simulations that suggest a higher value [23]. A real space 
renormalization group computation suggests that the upper critical dimension is infinite, in 
agreement with the numerically accurate (but not exact) formula for the critical exponents v\ 
controlling the width of the interface (i.e. tj = 2/(2D+ 1)). 

• The DLA (direct limited aggregation) problem. This problem is rather difficult and only a 
few attempts have been done to study the problem in large dimension There is the 

very interesting suggestion that the fractal dimension of the aggregate goes to D — 1 for large 
dimensions, but not much has been done in this direction. 

It is clear that each problem has its own peculiarities. Different approaches should be used. In 
all the cases, with the exclusion of the last one, the replica method can be used. In spin glasses in 
the low temperature phase the replica symmetry is spontaneously broken in the mean field approx- 
imation, while in the other cases it remains exact (with maybe the exception of electron glasses). 
In the next sections I will discuss in more details a few case. 

4. A difficult simple case. The random magnetic field Ising ferromagnet. 

In the nutshell, the random magnetic field Ising ferromagnet is a standard <p 4 theory with a 
quenched magnetic field. The probability of the field at given h(x) is 



with 



P((j)\h)^exp(-pS(<j)\h)) (4.1) 

S(<j>\h) = J dx(^(d(l>(x)) 2 + ^(t>(x) 2 + ^ 4 (x)-h(x)(l>(x) } . (4.2) 

Of course we must specify the probability distribution of the magnetic field h{x). We assume 
short-range correlations: 

P(h)=ex P (J-Jdxh 2 (x)^ , h(x)h(y)=g8(x-y) , (4.3) 

where the overline denotes the average over the h distribution. 

Going from positive X to large negative z we expect to find a phase transition: in the region of 
strongly negative % the expectation value of <p becomes different from zero and the Z2 symmetry of 
the model is spontaneously broken. A detailed and convincing analysis shows that critical behavior 
remains the same also in the limit /3 — > °°. A similar analysis show that near the critical point we 
have 



((j>(x)(i>(0))h~rix\h)r(0\h) . (4.4) 
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where (j)*(x\h) = (0(x))/, is the expectation value of the field <j>(x) in presence of a field h. 
Let us suppose that the solution of the following stochastic equation is unique 

= -A<j>(x) + Z(j)(x) + 0(x) 3 - h(x) = . (4.5) 

Let us call it (j)(x\h). In the limit j8 — s> °° = 0*(jc|/j). Under the previous assumptions it is 

possible to prove that then 

WJ¥W)d = (<K*)<M0))d-2 • (4.6) 
where {-)d-2 is the statistical expectation value with weight exp(— S[<j>\0]/g) (i.e. at zero magnetic 
field) in two dimensions less [10]. This dimensional reduction follows from an hidden supersym- 
metry of the stochastic differential equation. Indeed if the solution of the stochastic differential 
equation is unique, after some transformations one arrives to a field theory that is invariant under 
the 0(D\2) supergroup. From this result dimensional reduction follows. Unfortunately if the so- 
lution of the stochastic differential equation is not unique, one can prove the beautiful and useless 
identity: 

(0(jc)0(O))z>-2 = (0W0(O)> ss = £/„0„(jc|A)0 a (O|/i), (4.7) 

a 

where a label the solutions and I a is the Morse index of the (X th solution and the subscript SS 
denotes the supersymmetric results. 



If there are many solutions of eq. (4^5), Ea^a^ce (x\h)^) a (0\h) is very different from (0(jc)0(O))/j. 
Indeed we expect that solutions with lower free energy have an higher weight and we should have 



£„ I a exp(-5[0« |A] /g)4> a (x\h)<j> a (0\h) 



Unfortunately, when we approach the transition, z (i.e. the bare mass squared) becomes negative (it 
is always negative in the continuum limit), S((j) \h) is no more convex and we expect many solutions 
in the critical region. 

If we close our eyes on this problem, we can ask if dimensional reduction qualitatively works 
in three dimensions. The answer is no: the relation among the critical exponents of the random 
model in D dimensions and the non-random model in D — 2 dimensions, in particular 

Yrf{D) = Yp(D-2) , (4.9) 

cannot be satisfied, where Jp{D — 2) is the critical exponent of the pure system. The lower critical 
dimension is 2 for the random field problem (an exact result) and it is 1 for the non-random case. As 
far as 2 ^ 1 + 2, dimensional reduction fails badly. We should notice that the existence of multiple 
solutions is a phenomenon that cannot be seen in perturbation theory. Only non-perturbative effects 
destroy dimensional reduction and formally the 6 — £ expansion for random system coincides with 
the 4 — £ expansion for non-random system. A possible conjecture is 

7rf{D) = 7p(D-2) + 0(exp(-A/£)) . (4.10) 

Unfortunately it is not clear how to compute Jrf(D) — Jp(D — 2) . There are some sugges- 
tions (e.g. instantons), but no consensus. Moreover one could also take the view point that non- 
perturbative phenomena becomes important before we reach the transition and that the dimensional 
reduction formulae are not even approximatively valid |25]. The situation is quite confused and 
after nearly 30 years of work we are still far from a solution. 
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5. Anderson localization: Schrodinger equations for non-interacting electrons 

Lt us suppose that for some reasons we can neglect the interaction among electrons in a mate- 
rial. The properties of the electrons can be found by filling all the levels up to the Fermi level. At 
this end we have to find the solutions of the Schrodinger equations 

(-A + V(x))y(x)=Ey(x), U=-A + V, (5.1) 

where V is the potential due to the remaining part of the systems. The interesting case is when the 
potential V(x) is random (e.g. due the presence of impurities). Let us consider the simplest case: 
V(x) uncorrected from point to point with a Gaussian distribution: 

P(V)=expQ-J dxV 2 (x)^j . (5.2) 

Roughly speaking we can divide the states (i.e. the eigenvalues of H) into localized states 
and extended states. For localized states the spectrum is discrete and the wave function decays 
exponentially, e.g. as exp(— \x — xo\/%); the typical value of % is called the localization length. For 
extended states, the spectrum is continuum: in a finite volume the wave function of an extended 
state spread everywhere and it is order of 0(V~ l l 2 ) in each point of the lattice. 

It is possible to prove that for large negative E there are localized states. In dimensions D > 2 
there is an threshold energy El such that for E < El all states are localized, while for E > El all 
the states are extended. In dimensions D < 2 all states are localized. If the electrons fill the levels 
up to the Fermi energy, the conductivity is related to 

G 2 (x,E)=WMW, G(x|g) = (x|^|0)=£ ^^ (0) , (5.3) 

where E is complex (E = E + is): the real part (E) is the Fermi energy and the imaginary part (£) 
is proportional to the frequency at which the conductivity is measured. Indeed one finds that the 
conductivity a is given by 



a(e,E) oc £ 



2 dG(k,E) 



dk 2 



(5.4) 

k 2 =0 



where G(k,E) is the Fourier transform of G(x,E). In the following we will consider the limit 
£ — > 0. It is important to note that we cannot write the previous expressions for G(x\E) and Gi{x,E) 
directly in the limit £ = 0, because they would be undefined. Indeed 

\im {G{x\E + ie) - G{x\E -ie))=2niYy n {x)y n {0) 8 {E- E n ) . (5.5) 

This is similar to what is done in particle physics when we write G(p) = l/(p 2 — m 2 — ie). More- 
over in G2(x,E) there is a term that is proportional to 

I ^fSvfS " l -LWn(x)\ 2 W,m\ 2 8(E-E n ) . (5.6) 

It diverges (in distribution sense) when £ — > 0. The divergence disappears in the infinite volume 
limit when the states are extended, but it is still present in the localized phase. In the localized 
phase the zero frequency conductivity vanishes and the correlation functions do diverge. 



8 



Field theory for disordered systems 



Giorgio Parisi 



In the nutshell in the extended phase the correlation functions (e.g. G2 (*,£)) are long range, 
while in the localized phase the correlation functions are short range, but they diverge in the £ — > 
hmit. Essentially we have two possibilities for small e and k: 

The first is realized in the extended phase and the second is realized in the localized phase. 

We would like to understand the behavior of the correlation length and of the conductivity 
near the localization transition, i.e. near El. At this end it is expedient to use a field theory 
representation. We introduce n real fields a and n real fields <p a . Eventually n — > 0, like in quenched 
QCD. The theory in the e — > limit has an Euclidean action invariant under the non-compact 
group 0(n,n). This group is defined as the set of transformation that leaves invariant the following 
quantity: 

QM^ I (fa~~€) • (5-8) 

a=l,n 

The explicit form of the action is 

y</[0]exp(-S[0]), (5.9) 

where S[<f>] is given by: 

ifd d xt £ ((d,U^) 2 -(d^ a (x)) 2 -is^ a +^)+EQ 2 (x) + g{Q 2 (x) 2 )) ) (5.10) 

The factor i in the action is crucial: the term £ (d)% + 0^ explicitly breaks the 0(n,n) symmetry, 
but is presence is necessary in order to make the functional integral convergent. 

What happens to the 0(n,n) symmetry when e — )■ 0? This symmetry is broken in the non- 
localized phase and in some sense it becomes an exact symmetry in the localized phase. The 



peculiar behavior shown in eq. (5.7) is characteristic of non-compact symmetry groups: the Green 
functions are infinite in the unbroken symmetry phase. In this representation localization is a non- 
perturbative phenomenon (it arises as the effect of instantons). In order to compute the properties 
of the theory, we should study the critical behavior of a strongly interacting multi-instantons gas 
and this is not an easy task. Our inability to control the localization transition in this representation 
implies that the upper critical dimension in not yet known (may be it is infinite!). At low dimensions 
things become simpler. We have seen that there are no extended states in D = 2 (it is a variation of 
a well known no-go theorem). There is an other field theory representation (a non-linear a-model), 
where a transition is present in perturbation theory from the beginning and the exponents are know 
in 2 + £ dimensions. In conclusion electron localization is an example of a theory where the lower 
critical dimension is know, but nothing sure is known about the upper critical dimensions. 

6. Spin glasses 

Many progresses have recently done in spin glasses: theory, experiments, simulations and even 
theorems! The simplest Ising spin glasses has the following Hamiltonian: 

H j{ a ) = \ L JikOiOk (6-1) 

1 ik=l,N 
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where the a are Ising spins (i.e. ±1), the i and k are neighbors. The J's are random (e.g. ±1) and 
they are defined on the links of the lattice. If a magnetic field is present we must add to the Hamil- 
tonian a term proportional to hO\. Let us present an alternative definition (in the continuum) that 
may be more appealing: they expose some of the symmetries of the problem. Let us consider a 
gauge group G, a gauge field A^(x) and a gauge transform g(x). In the Landau gauge we want to 
minimize 

J dxlr(A 8 )1 {x)) 1 =H A {g). (6.2) 
We can consider (for a given gauge field) the probability distribution of g given by 

exp(-j3tf A (g)) . (6.3) 

This is a generalization of spin glasses, where g and A correspond to a and /. We obtain the 
previous definition when the gauge group G is Z2, we are on the lattice and we consider the strong 
coupling limit. In many cases Gribov ambiguity tells us that Ha (g) has many minima, therefore 
we should not be surprised to known that Hj(a) has an exponentially large number of minima. We 
can ask many questions. Which is the form of fl^(g)? How different are the minima? Which are 
the properties of set of minima (e.g. of their mutual distance?) Briefly which are the properties of 
the (free) energy landscape? 

This problem can be studied in the mean field approach where one can obtain exact results [^] 



(partially corroborated by theorems [|26|]). It turns out that there is low temperature phase where 
the system has an exponentially large number of equilibrium states that are local minima of the free 
energy. Each state has a very large (exponentially large) mean life and at equilibrium the system 
travels from one state to an other state. 

The system has a very slow approach to equilibrium. There is a fast dynamics inside the 
states and a very slow dynamics due to the barriers that separate the states (like punctuated equi- 
libria in biology or in the landscape theory for strings). A peculiarity of this dynamics is the 
presence of violations of the fluctuation dissipation relations in slightly off-equilibrium dynamics. 
Numerical simulations give ample evidence of these new fluctuation-dissipation relations. These 
new fluctuation-dissipation relations have been observed in some difficult experiments (notably by 
Herisson and Ocio in spin glasses [§]). 

May be the most characteristic phenomenon is the presence of two different magnetic suscep- 
tibilities that correspond to different paths in the temperature magnetic field plane. In the nutshell 
we have: 

• The magnetic susceptibility that correspond to the response in one state (Xlr) that is mea- 
sured by adding an infinitesimal magnetic field at fixed temperature. 

• The true equilibrium susceptibility, where one includes the effects of changing the equilib- 
rium state when we change the magnetic field. It can be experimentally approximated by the 
field cooled magnetization Xfc- 

The naive predictions of mean field theory and the experimental data agree in an impressive manner. 

What happens in finite dimensions (we are particularly interested to the three dimensional 
case)? We would like to characterize the infrared (large scale) behavior in the low temperature 
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phase and find in which dimensions the low temperature phase may exist. Numerical simulations 
are in good and detailed agreement with the predictions of a generalized replica approach in D > 3 
(at least in zero magnentic field), while it is clear that in D = 2 there is no transition (and the lower 
critical dimension is greater than 2). 

By interpolating various exponent as function of the dimensions one obtains the estimate Dl = 
2.493(7) Jl9|]. An explicit theoretical argument based on the evaluation of the interface energy 
gives the unexpected value Du = 2.5. Therefore we have a very good agreement among the theory 
and numerical simulations. However the value Dl = 2.5 is strange. In the usual computation the 
lower critical dimension can be related to the infrared divergence of some simple diagrams. In 
scalar theories the spontaneous breaking of the 0(N) symmetry implies the presence Goldstone 
Bosons (the corresponding propagator being 1 /k 2 ). The interaction is proportional to k 2 (soft pions 
theorems) and infrared divergences pile up in dimensions less or equal to 2. In spin glass we are 
unable, as we shall see later, to do such a simple computation. It is difficult to think which simple 
diagrams could diverge starting from 2.5 dimensions. We would like to answer to the following 
questions. Can we make a similar simple perturbative argument in spin glasses? Which is the low 
k behavior of the propagators in D=3? We may get it from numerical simulations, but which are 
the theoretical predictions? 

Mean field theory was solved by replica approach. Before presenting the field theory con- 
structed as an expansion near the mean field theory we must spend a few words in describing 
some of the unusual characteristic of the mean field. As usual we need to compute the average 
value of the logarithm of the partition function (i.e. the connected vacuum to vacuum diagrams): 



F = -ln(Zj), (6.4) 

where Zj is the J dependent partition function and the overline is the average over J. It turns out 
that it is simpler to compute first 

log(Zf) 

F n = • (6.5) 

n 

At the end we recover the value of F using the formula 

F = limF„. (6.6) 

This approach leads to a strange mathematics: one introduces a is n x n matrix Q a b (a,b = l,...n) 
and the symmetry group is S n (the permutation group of n elements); eventually n goes to 0. In this 
way we construct a free energy F [Q] and we obtain 



dF 

F = F[Q*\ , where 



dQa,i 



= 0. (6.7) 

Q=Q* 



As a first approximation we suppose that the group of symmetry of the system (S n is the 
permutation groups of the n replicas) is spontaneously broken to it subgroup S n / m ® {S m ) n l m : the 
replica symmetry is spontaneously broken (it can be seen that spontaneous breaking of replica 
symmetry correspond to the existence of an infinite number of equilibrium states). When n — > 0, 
Sn/m ~~ ^ So so that So is a subset of So and the unbroken subgroup contains So. This fact implies 
the possibility of breaking the permutation group in an infinite hierarchical way. The permutation 
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group So is promoted to a continuous group. Indeed one can introduce infinitesimal permutation 



and write the Ward identities associated [jig, |17|]. At the end of the day the matrix Q is parametrized 
in terms of a function q(z) defined on the interval 0—1. One finally finds a free energy F[q] (it a a 
functional of q(z)) that must be maximized (not minimized!) []|, 27]. Talagrand theorem [ |2"6| ] tells 



us that this baroque construction gives the exact result in the Sherrington Kirkpatrick model. 

The perturbation around mean field theory can be formally constructed. There are no diffi- 
culties in the high temperature phase (the unbroken symmetry phase is conventional). We know 
the upper critical dimension (i.e. Dy = 6) and we can control the critical exponents in the 6 — £ 
expansion (three loops have been computed). The problems come in the low temperature phase 
where we do not have a good control of the infrared behavior. The perturbative expansion has been 
constructed, but the evaluation of the results is difficult. As usually one starts by computing the 
propagator at tree level: 

G a , b;c ,d(x) = (Q a ,b(x)Qc,d(Q))c , (6-8) 

(we will denote by G a j,- C ^{k) the Fourier transform of G a} b; C ,d(x). 

The propagator depends on four indices and in the replica broken phase it can parametrized in 
terms of functions of three variables zi,ZiiZ3. The final form of the propagator can be obtained by 
solving six coupled integral equations [ |l6| , 17]. Depending on the choice of the indices (a,b,c,d) 



G a ,b;c,d (k) has a behavior near k = of the type 

0(1), 0(l/k 2 ), 0(l/k 3 ), 0(\/k 4 ) . (6.9) 

These long-range correlations (i.e. massless particles) are related to the existence of Goldstone 
modes of a symmetry that is spontaneously broken in the low temperature phase and some of 
them are consequences of the Ward identities [[if], 17] associated to the infinitesimal permutations. 



Moreover there is an infinite number of poles at k 2 = —mf, where the index i goes up to infinity 



and the m,-'s accumulate toward zero Q16Q. The one loop corrections can be written in an explicit 
way. If we proceed in a naive way we find at the end 0(1O 4 )- 0(1O 6 ) different integrals and the 
computation is nearly unfeasible. In the infrared region there are very strong cancellation among 
the various terms (e.g. [|28[]) so that a very careful evaluation is needed. It clear that one should 
regroup the various terms in a smart way, may using Ward identities. We have not been sufficiently 
smart. 

The final effect of this complicated behavior is that we are not able do construct the appropriate 
non-linear sigma model that should be used near the lower critical dimensions. The Ward identities 
of the symmetry have been identified and they have been shown to be the responsible of many 
important cancellations, but we are still far from writing the final effective theory [|i"6|]. Maybe we 
are very near, but we miss at least one crucial ingredient and we are stuck. 



7. The topological expansion: Fat Diagrams 



In this section I would only mention a new idea that has bee put forward by Efetov [ 29, 30, 3 1 ] 
about twenty years ago. It looks very promising; however not too much progress has ben done 
because the needed computations are feasible, but technical difficult. 
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This approach can be defined only for a theory on the lattice. The typical Hamiltonian is of 
the form: 

i,k 

where the sum is done over the nearest neighbor points, the Hamiltonian may be non-translational 
invariant (has happens in a random system) and the field 0,- may be a vector with many components. 
Our aim is to find out the relevant degrees of freedom in the continuum limit, i.e. near the phase 
transition where the the correlation length becomes infinite. There are many ingredients that enter 
in the approach: a) the Bethe version of mean field theory that becomes exact on the Bethe lattices; 
b) the control of the one dimensional version of the model; c) the high temperature expansion in 
lattice theories. 

Let me try to summarize the basic ideas. In mean field theory, one neglect the loops and 
one looks for a self-consistent solution of the equation of motion. This corresponds to neglect 
correlations and to impose self-consistent boundary conditions. Mean field theory becomes exact 
in the case where the lattice is locally loopless (i.e. it has the topology of a tree): loops are present, 
but they have a length that goes to infinity with the volume (this happens on Bethe lattices). The 
presence of large loops accounts for the need of imposing self consistent conditions for the vacuum 
expectation values. 

In the topological expansion one starts from the mean field theory, where the lattice is locally 
a tree. In order to describe the approach it is convenient to recall the particle representation in field 
theory. In D dimensions the free propagator (i.e. l/(k + m ) in momentum space) can be written 
in configurations space as: 

G(x) = J ds r yV(x,s)exp(-m 2 s) , (7.2) 

where JV (x,s) °< s~ D / 2 exp(— x 2 /s) is the probability of a path of length (internal time) s going 
from to x. This representation is well known in the continuum, but it can be extended also to the 
lattice. 

A one loop contribution for the propagator in a 3 theory can be written in a similar way as 

J d D yd D z.G(x-y)G(x-z) 2 G(z) = J d D yd D z J dsids 2 ds 3 ds A (7.3) 
exp(-ra 2 (si + s 2 + S3 + Si))J/ '{x-y,s\)jY '(y - z„s 2 ),yV '{y - z„s 3 ),yV \z,s 4 ) ■ 

In the fat diagram approach we will derive similar formulas in a different context: we have 
to consider a theory on the lattice in one dimensions in a fluctuating background and at the end, 
in eq. (fO|), we will substitute exp(— m 2 s) with the correlation function on an one-dimensional 
model. We substitute in eq. (7.3) exp(— m 2 (s\ + s 2 + sj, +S4)) with the correlation function on an 
one-dimensional model, with the topology of the corresponding Feynmann graph. 

Let us see how to proceed in details. If we neglect loops we would naively write the following 
equation for the probability distribution of the field : 

II (/ d(j> k P((l> k ) e xp{-H((j>,^ = (| JfP(f)exp(-//(ff )) , (7-4) 

where z in the coordination of the lattice (e.g. z = 2D in an hypercubic lattice). However the 
previous formula is not correct because the fields in two nearest neighbour points are correlated. 
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The correct results for a loopless system is given by the Bethe approximation. We call B(<j>) the 
probability distribution the field with z — 1 neighbours (we remove one spin and we consider the 
probabability distribution of the spins around this cavity). After some work get the equations: 

B(0)oc (yV*(f)exp(-ff(M')) Z > (/W)exp(-ff(ff)j Z . (7.5) 

Let us compute the correlations in this framework. Let us suppose that there is an unique path 
from point x to y of lentgh L. We can compute the probability 

^ L) (0(x),0(j)) =P(0(*),0(y))-P(0(*))P(0(y)) • (7-6) 

by solving the one dimensional problem where each site interacts with z — 2 sites outside the unique 
one dimensional path. Let us assume that P c (<j>(x),(j>(y)) is small, so that if there are many path 
they do not interfere. We finally get the zero loop contribution: 

P°(«K*),<K)0) =£^(L,x-j)P c (L) (0(x),0(j)) . (7.7) 

L 

The effects of the loops can be included in a systematic way, by considering the contribution 
of lattice region with high genus. In this way one obtains new objects (that we can call them fat 
diagrams): they look like usual diagrams, but they contain the resummation of non-perturbative 
effects in one dimensions. Their computation is rather complex, but it seems to be feasible. At the 
end of the day we get a genus expansion: 

P c (0(*),*(y)) = £ P*(<!>(x),<!>(y)), (7.8) 

g=0.~ 

where P§ (0 (x) , (j) (y) can be computed on a lattice of genus g. 
8. Conclusions 

We have seen that there are many simple disorder systems where we do not have a good control 
of physics (e.g. we do know the behavior near the upper critical dimension, if any, and near the 
lower critical dimension): 

• Equilibrium systems, e.g. ferromagnets with a random field, localized electrons, directed 
polymers, spin glasses (both with and without a field) . 

• Non-equilibrium systems, e.g. surface growth (KPZ equation), diffusion limit aggregation 
(DLA). 

In particular we do not control the large scale behavior of spin glasses in the broken symme- 
try phase. Many of these problems stand for twenty years or more. We badly need new ideas, 
new tools, techniques and also smart people trying to solve them. Fat diagrams is a promising 
possibility. I hope that they will fulfill the promises! 
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